Voter models with conserved dynamics 
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We propose a modified voter model with locally conserved magnetization and investigate its phase 
ordering dynamics in two dimensions in numerical simulations. Imposing a local constraint on the 
dynamics has the surprising effect of speeding up the phase ordering process. The system is shown 
to exhibit a scaling regime characterized by algebraic domain growth, at odds with the logarithmic 
coarsening of the standard voter model. A phenomenological approach based on cluster diffusion 
and similar to Smoluchowski ripening correctly predicts the observed scaling regime. Our analysis 
exposes unexpected complexity in the phase ordering dynamics without thermodynamic potential. 
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How do systems of interacting particles order, and 
what are the mechanisms by which microscopic in- 
teraction generates macroscopic structures? These 
questions have first been asked in biology [T] and 
then in condensed matter physics, magnetism and 
statistical mechanics, where kinetic ordering dynam- 
ics has received considerable attention [21 EH- Or- 
dering processes have more recently also become of 
interest in interacting agent models of social pro- 
cesses [H [5]. Answering our opening questions re- 
quires an off-equilibrium approach starting from the 
microscopic dynamics itself - a challenge of consid- 
erable intricacy. Our understanding of processes out 
of equilibrium is still limited, and developing a more 
complete picture is a focus of current research. 

Processes of kinetic ordering are traditionally ad- 
dressed in the context of simple lattice models. In 
condensed matter physics their dynamics is usually 
based on the minimisation of a thermodynamic po- 
tential, coupled to a heat bath and obeying detailed 
balance. Such models allow for a characterization of 
different types of ordering. There is a distinction be- 
tween models with and without a local conservation 
law. Glauber dynamics and the kinetic Ising model 
of magnetic systems have no conservation law, their 
characteristic length scale grows algebraically with 
time, L{t) ~ f™, where the coarsening exponent is 
found to take the value n = 1/2 in phenomenological 
theories [H[3]. In the conserved Kawasaki dynamics, 
on the other hand, one finds n = 1/3 [21 IS]- This 
type of dynamics is used to describe alloys or binary 
liquids. In both models the coarsening dynamics is 
driven by surface tension. 
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A second class of dynamics is defined by systems lack- 
ing a thermodynamic potential or energy function. 
Such models are frequently motivated by human so- 
cial interaction ,5J. This type of model, frequently 
based on pair-wise interaction, has also been used to 
study autocatalytic chemical reactions [5] and bac- 
terial populations [7 . In lattice models with Z2- 
symmetry and absorbing states domain coarsening 
can proceed in the absence of surface tension, driven 
by interfacial noise [HI IS]. The main representative 
of this class is the voter model (VM) [TU], in which 
spin-like variables align with the state of a randomly 
chosen nearest neighbor. This dynamic does not per- 
mit a locally conserved quantity, even though global 
magnetization is conserved in the ensemble. The 
coarsening dynamics of the VM is algebraic (with 
n = 1/2) in spatial dimension d = 1, and logarith- 
mic in d = d c = 2, that is, L(t) cx 1/lnt QH[T2|. For 
d > d c the interfacial noise becomes irrelevant and 
the infinite system does not order. 

The purpose of our work is to investigate the phase 
separation dynamics in models of the voter type, 
but with local conservation of the order parameter. 
Specifically we study spin-exchange processes based 
on interactions of pairs of particles. Similar to the 
celebrated Schclling model of segregation [3] such dy- 
namics can be motivated by social processes. This is 
not the main objective we wish to pursue though. We 
aim to systematically work towards a more complete 
picture of the possible types of ordering dynamics in 
off-equilibrium particle models. Specifically we in- 
troduce variants of the VM with local conservation 
laws. While some VM modified along these lines do 
not order at all, others show an effective algebraic 
domain growth. We provide a phenomenological un- 
derstanding of this behaviour invoking a combination 
of two different mechanisms: (1) a linear instability 
triggering the formation of compact faceted patterns, 
and (2) a process similar to Smoluchowski ripening 
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[15] , responsible for domain coarsening in the long- 
time regime. The latter occurs through cluster co- 
alescence driven by surface diffusion and follows an 
algebraic law L(t) ~ i 1 / 5 . Interestingly, imposing the 
local constraint does not slow down the dynamics, 
but speeds it up relative to the logarithmic ordering 
of the standard unconstrained VM [14] . 
Microscopic Dynamics. Consider binary variables 
Si G {+1,-1} defined on sites i of a d-dimensional 
lattice. At each time step two neighboring sites i and 
j are chosen at random. If Sj ^ Sj then the two spins 
are exchanged with probability 
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neighborhood of i. These dynamics can be under- 
stood as follows: each of the two sites i,j randomly 
selects one of their remaining neighbors, k G di\ j 
and I £ dj \ i. If sg = and Sk — Sj, then Si and 
Sj are exchanged. A more detailed motivation can be 
found in the accompanying EPAPS document [15], as 
well as a discussion of other variants of the model. 
We consider unbiased random initial conditions in 
d = 2 dimensions. The temporal behavior of the 
density of interfaces (neighboring pairs of opposing 
spins) p(t), shown in Fig. [ija), reveals a very slow 
initial dynamical regime, in which clusters of spins 
with same orientation are formed. At long times 
the system enters a regime of algebraic coarsening, 
in which the density of interfaces decays according to 
p(t) ~ t~ x . This behaviour is similar to that seen 
in curvature-driven phase separation, and in contrast 
to the usual VM in d — 2 which exhibits logarithmic 
coarsening. In this dynamic scaling regime the sys- 
tem has a single characteristic length scale, L(t), and 
the order-parameter correlation function assumes a 
scaling form C(r,t) — f(r/L(t)). The typical length 
scale can be extracted from the first zero of the cor- 
relation function, C (r, t) and is found to grow alge- 
braically, L(t) (x t n , see Fig. [ijb). From our data the 
coarsening exponent is found as x ~ n ~ 0.21, and 
therefore different from standard the exponent 1/3 
in curvature-driven coarsening with locally conserved 
order-parameter [2J. Snapshots from simulations of 
the conserved VM are shown in Fig. IlJ(e), with time 
the system develops faceted domains with straight 
and elongated diagonal boundaries. 
Pattern formation and mean-field approach. Unlike 
standard VM the stochastic dynamics of the present 
model cannot be solved exactly due to the non- 
linearity of the transition rates. Following classical 
approaches [3j, we neglect correlations between dif- 
ferent sites and derive mean-field evolution equations 
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the product of their averages (sj) • • ■ (sj) . We will use 
the variable ifi to identify the mean-field approxima- 
tion of the average spin value (s,). In a continuous- 



ly 



! v 


■ 


1 


■ 



10 10' 10 10 

t 




t = 10 5 t = 10 6 t = 10 7 



FIG. 1: Voter model with microscopically conserved dy- 
namics. Panels (a) and (b) show the time-evolution of 
the density of interfaces and correlation length respec- 
tively, panels (c) - (e) snapshots of the kinetic ordering 
process. Dashed lines in (a) and (b) indicate algebraic 
growth laws with exponent 0.21. Data is from simula- 
tions of a 500 x 500 system. 



time limit, the evolution equation for the local mag- 
netization (pi becomes 
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resulting in a lengthy expression on the right-hand- 
side, reported in detail in the accompanying EPAPS 
document [TS]. These mean-field equations satisfy 
the conservation law Lpi = 0. The mean-field de- 
scription can be used to explain the process of pattern 
formation. Linearizing ^ around the homogeneous 
solution ipi = Vz, one finds 



(3) 



where V 2 ^^ is the (lattice) Laplacian. The present 
model is therefore characterized by a linear instability 
of type II in the Cross-Hohenberg classification [T6] 
with a band of unstable wave numbers < k < 1/ Vd 
with maximum growth rate at fc max = l/y2d. A lin- 
ear instability generated by terms of the type as in 
(J3j) is typical of coarsening with conserved dynamics, 
described by the Cahn-Hilliard equation (2J. These 
terms have a different origin in our model though. 
In the Cahn-Hilliard equation the term proportional 
to — V 2 (/? comes from linearizing the thermodynamic 
potential; in the present model it is instead generated 
dynamically by exclusion of spin j from the sampling 
of i's neighborhood and i from the sampling of j's 
neighborhood. This restriction prevents bulk diffu- 
sion, i.e. the migration of isolated defects, which 
is present in the Kawasaki spin-exchange dynamics 
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FIG. 3: (Color online) Periphery diffusion: The exchange 
process of the two spins marked by the red (grey) frame, 
along with their respective neighbors (marked by dots), 
leads to a rearrangement of kinks at the interface between 
two clusters. This may result in a diffusion dynamics of 
the centre of mass of clusters (cluster diffusion). 



FIG. 2: (Color online). Interface density from a numerical 
integration of the mean-field dynamics, Eq. pj, and from 
the microscopic model with f2 spins per site. Data is from 
single simulation runs of a system of size 200 x 200. 



even at T = [3] . The influence of the restriction be- 
comes progressively less important in higher dimen- 
sions, and the linear instability is expected to disap- 
pear for d — y co. If we remove the site restriction, 
only the surface diffusion term, — V 4 <£>, survives in 
Eq. ([3]) , and phase separation does not occur [T5] . 

Although the mean-field approach describes the pat- 
tern formation process qualitatively to some extent, 
it does not capture the correct long-time coarsening 
regime. This is not surprising, in the example of the 
standard VM the mean-field dynamics is simply given 
by diffusion equation, and it does not capture the 
essential features of the dynamics in low spatial di- 
mensions (e.g. d = 2), where the fluctuations due to 
multiplicative interfacial noise are relevant. The con- 
served VM we discuss here is more complicated, nev- 
ertheless the differences between the dynamics of the 
original microscopic model and the mean-field equa- 
tions can again be traced back to the discreteness of 
the degrees of freedom, i.e. to multiplicative noise. 
We have considered a modified model in which each 
lattice site contains VL spins, similar to what was pro- 
posed in [17) . In a given update step a spin interacts 
with randomly chosen spins in the neighboring sites, 
following the microscopic rules of the conserved VM. 
For f2 — > 00, the dynamics is deterministic and given 
by the above mean-field equations. The data in Fig. [2] 
show that for f2 large but finite, the temporal behav- 
ior is very similar to the predictions of the mean-field 
equations in an initial brief coarsening regime gov- 
erned by surface diffusion (also referred to as 'bidiffu- 
sion' [3]) and with algebraic coarsening with n pa 1/4. 
Then the pattern formation process takes place and 
the dynamics finally reaches a structured configura- 
tion of patterns, with only a slow further coarsening 
process or potentially none. Considering the case of 
small ft ~ 1, instead, the system develops an alge- 
braic coarsening regime with exponent n ~ 0.2, which 



persists in the long-time regime. 

Coarsening by Cluster Diffusion. At this point it is 
clear that the observed asymptotic coarsening regime 
is an intrinsic effect of individual discreteness, and 
that it cannot be captured by simple mean-field equa- 
tions. Instead we follow a phenomenological ap- 
proach, relating to a phenomenon known in material 
science as Smoluchowski ripening [13] . to obtain fur- 
ther insight into the process of cluster aggregation. 
In absence of bulk diffusion, we can exclude Ostwald 
ripening [18] . therefore we expect a slow coarsening 
process, with n < 1/3. Horizontal and vertical do- 
main walls are blocked, while other interfaces can still 
move by means of a process known as 'periphery dif- 
fusion', that is a slow diffusional movement of kinks 
on the boundary of a domain (see the illustration in 
Fig. [3| . This movement causes a slow effective dis- 
placement of the center of mass of the clusters, i.e. 
an effective cluster diffusion process. Larger clusters 
move more slowly because the effective diffusion co- 
efficient of the cluster depends on the ratio between 
the surface and the volume of the clusters. The form 
of the resulting diffusion coefficient of clusters has 
been derived in the literature, see e.g. [TS]- In or- 
der to move the center of mass of a cluster by a dis- 
tance S c in two dimensions, we need order L 2 spins 
to move by the same distance. This corresponds to a 
single spin moving over a distance S p oc L 2 5 c . On the 
other hand, the rate with which moves of the centre 
of mass of a cluster occur is proportional to that of 
diffusion events on its surface, i.e. T c oc T p L. Given 
that the single spin diffusion coefficient is D p oc S p T p , 
we obtain that the diffusion coefficient of a cluster 
of typical scale L decays as D C (L) oc L~ 3 . In or- 
der to find the coarsening exponent we next consider 
the temporal evolution of the average cluster size by 
means of the Smoluchowski equation for cluster coag- 
ulation processes. This mean-field approach is usu- 
ally considered qualitatively correct also in two di- 
mensions, Kang and Redner [20] have shown that the 
upper critical dimension is d c — 2. We therefore ex- 
pect logarithmic corrections at most. We follow the 
derivation proposed by Kandel [31], and introduce 
the size (area) of a cluster in d dimensions, s. The 
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FIG. 4: (Color online). Scaling behavior of the cluster- 
size distribution P(s,t). Main panel shows the collapse 
plot obtained using the predicted scaling relation Eq. Q , 
the inset shows the raw data for P(s,t). Data is from 
10 runs of a system of size 400 x 400 (running average 
performed for data in main panel, see also [150. 

time-dependent density of clusters of size s per lat- 
tice site is denoted by P(s, t) [22 . One then finds the 
following scaling relation in the long-time limit |21j 

P(s,t) = t~ a f [s/sav(*)] , (4) 

where the average cluster size s av (t) = 

£ s s^fs, i)/ scales as s av(i) ~ tP. The 
scaling exponents are predicted to be a = 2/(£ + 1) 
and (5 = + 1), where C characterizes the scaling 
of the diffusion constant of clusters, with the cluster 
size, D c (s) ~ s~ c 21J. The scaling of Eq. Q is 
confirmed by simulations of the microscopic model, 
see Fig. |4j adding weight to the hypothesis that 
a process similar to Smoluchowski ripening is at 
work in our system. In the present case we have 
D C {L) ~ L -3 , i.e. £ = 3/2. Therefore we expect that 
the asymptotic behavior of the microscopic model 
is characterized by a coarsening regime s av ~ t 2 ^ 5 . 
Given that L is the unique length scale the hy- 
pothesis of Smoluchowski ripening leads to L ~ t n 
with n = fi/d = 1/5. The data shown in Fig. |T| is 
consistent with this phenomenological prediction. 



We note though that the facetted configurations 
shown in Fig. [I] differ from those seen in physical 
systems undergoing Smoluchowski ripening |23j . In 
these systems one species of particles is frequently 
more abundant than the other, which may explain 
the difference in the resulting structures [24] . 

Globally conserved dynamics. To determine if the 
modified scaling behaviour we observe is genuinely 
an effect of local (microscopic) conservation, as op- 
posed to conservation per se, we considered a version 
of the dynamics with conservation only at a global 
level [T5]. We find logarithmic coarsening as in the 
standard VM [T5] , in-line with earlier work indicating 
that conservation laws at long-range do not affect the 
coarsening dynamics of non-conserved models [25] . 
Conclusions. We have put forward a voter-like 
model with strict microscopic conservation of mag- 
netization. The spin-exchange rule introduces non- 
linearities that turn the logarithmic growth of the 
standard VM into an algebraic coarsening process, 
driven by diffusion at the surfaces of clusters and by 
coalescence of clusters, similar to what is known as 
Smoluchowski ripening. Simulation results are con- 
sistent with a growth law L ~ i 1 / 5 , as predicted by a 
phenomenological theory. These findings extend the 
complexity and variety of the ordering of interact- 
ing particle models out of equilibrium. In particular 
our analysis exposes unexpected behavior in models 
with conserved dynamics lacking a thermodynamic 
potential. This, we think, can be relevant not only 
for non-equilibrium statistical physics and biological 
processes, but also for models of segregation in the 
social sciences. 
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